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ABSTRACT 



We report the characterization of bright, compact features in the cosmic 
microwave background radiation (CMBR) detected during the June 1992 
and June 1994 balloon flights of the Medium Scale Anisotropy Measurement 
(MSAMl-92 and MSAMl-94, respectively). Spectral flux densities are 
determined for each feature at 5.7, 9.3, and 16.5 cm~^. No viable counterparts 
for these features were found in source catalogs at 5 GHz or at 100 fim. The 
measured spectrum of each feature is consistent with a temperature fluctuation 
in the CMBR. The existence of these features is consistent with adiabatic 
fluctuation models of anisotropy in the CMBR. 



Subject headings: cosmic microwave background — cosmology: observations 



1. Introduction 

Recent observations of the cosmic microwave background radiation (CMBR) have 
found statistically significant anisotropy at 0?5 angular scales (see [White, Scott, and Silk 
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1994, and references therein). In earlier Letters ( Uiieng et at. 199i , hereafter Paper I, and 
Uheng et al. 199"6| , hereafter Paper II), we reported on observations of CMBR anisotropy 
that included two particularly bright features. The presence of bright, unresolved sources 
suggested the possibility of foreground point-source contamination, or of non-Gaussian 
fluctuations in the CMBR. In this Letter, we present our analysis of these features and a 
discussion of other work studying their properties. 



2. Instrument and Observations 

The MSAMl instrument has been described in detail elsewhere (|I''ixsen et at. 19^BD; 



only the most relevant features are given here. MSAMl is a balloon-borne off-axis 
Cassegrain telescope with a nutating secondary mirror and a four-band bolometric 
radiometer. The detectors have center frequencies of 5.7, 9.3, 16.5, and 22.6 cm~^ each with 
~ 1.5 cm~^ bandwidth. The beam size of the telescope is 28' FWHM, roughly independent 
of frequency. The secondary executes a horizontal four-position square- wave chop (i.e., 
center, left, center, right) with an amplitude of ±40' on the sky. The telescope is pointed 
using a combination of torque motors and reaction wheels, referenced to a gyroscope. The 
absolute sky location is determined by star camera observation. 

The data discussed here were collected during two flights of the MSAM instrument 
from Palestine, Texas; the first on 1992 June 5 (MSAMl-92), the second on 1994 June 2 
(MSAMl-94). These observations have been described in Papers I and II, and are only 
summarized here. In each flight, the antenna patterns for all four bands were measured 
by rastering over Jupiter. We determine the location of the antenna pattern within the 
star camera's field-of-view to within 2' (post-flight) by simultaneously observing Jupiter in 
the star camera and the main telescope. Gyroscope drift increases our absolute pointing 
uncertainty to ±2f5. We calibrate our spectral flux density measurements using the 
measured surface temperatures ([Griffin et al. 1986|) together with Jupiter's apparent 
diameter at the epoch of flight ( [Hagen and Boksenberg 1991| , [Donat III and Boksenber^ 



1993|) . The calibration has a 10% overall uncertainty, with relative errors of 1-4% between 



bands. 

The CMBR observations were made at roughly constant declination 5 ~ 82°, with the 
central beam covering right ascensions litA-lQ^Q and 17^^2-20^3 during the MSAMl-92 
flight. This corresponds to ~ 6.6 square degrees of sky coverage. The goal of the MSAMl-94 
flight was to observe the identical field as MSAMl-92; this was nearly achieved, with an 
average declination difference of ~ 10'. The coverage in right ascension was somewhat less, 
extending from 15'?0-17':^0, and 17^^3-20^0. 
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Except for §^ all results reported in this Letter are based on the MSAMl-92 
observations. 



3. Data Analysis 

The initial data reduction used in this analysis has been reported in Papers I and II; 
only the salient features are reviewed here. Note that this Letter is based on the reanalyzed 
version of the MSAMl-92 data including the full instrument covariance matrix (Inman 



et at. 19961 ). Two nearly statistically independent demodulations are formed from the 
data. The single-difference demodulation subtracts the right beam from the left beam, and 
ignores the central beam data; the double-difference demodulation subtracts the average 
of the left-t-right beams from the central beam. The single-difference demodulation has an 
antisymmetric two-lobe response pattern, while the double-difference demodulation has a 
symmetric three-lobe pattern. The data are binned by position and angular orientation 
of the demodulated antenna pattern on the sky. To remove instrumental drift, a slowly 
varying function of time is fit simultaneously with the bins. The bin size is 0?12 in position, 
and 10° in angular orientation. 

The dominant sources of brightness contrast expected in our bands are anisotropy 
in the CMBR, and thermal emission from warm (~ 20K) interstellar dust ( [Weiss 1980| , 



Bennett et al. 1992 ). Before searching for compact features, we use one of our spectral 
degrees of freedom to make a simple correction for the presence of dust emission. Unlike 
our previous analyses for CMBR anisotropy, here we do not impose a full spectral model 
on the data. Rather, after a minimal correction for dust contamination, we model the 
spatial structure of several bright, compact features, and study the spectral nature of those 
features. 



3.1. Foreground Dust Subtraction 

We estimate and subtract the contribution of foreground interstellar dust emission 
using the 22.6 cm~^ band, assuming a z/^'^ emissivity Rayleigh- Jeans spectrum (the noise 
estimates in the remaining three bands are adjusted to reflect this subtraction). In making 
this correction, we assume that the 22.6 cm~^ band is an uncontaminated monitor of dust. 
Varying the exponent of the power law emissivity of the dust by ±0.2 has a negligible effect 
on the final results, compared with the statistical uncertainties. The dust correction is most 
important for the 16.5 cm~^ band, and has little effect at the lower frequencies. 
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3.2. Search for Unresolved Features 

After removing the estimated dust contamination, we search the remaining 3 frequency 
bands for unresolved features using an F-statistic to measure significance. Each band is 
searched independently, and then the source lists are compared. Since the telescope is 
scanned at roughly constant declination, we parameterize the location of the features as 
a function only of right ascension, and constrain the central declination of the feature to 
follow a slow linear function of right ascension. 

We begin by calculating the raw of the data in each band, 

where tj are the binned data with associated noise covariance Vij, and the sum is over both 
demodulations as well as bins. Oj is one of four arbitrary offset parameters corresponding 
to the two demodulations and two portions of the flight, appropriate to datum i. Then, as 
a function of right ascension a, we flnd the flux density Sp{a) that minimizes 

= Y.[t^-0,- S,{a)B,{a)] [t, - O, - S,{a)B,{a)] , 

where Bi{a) is the value of the beam map (single or double difference, as appropriate) 
centered on observation point z, evaluated at right ascension a. The signiflcance of the 
fitted fiux density is determined by forming the statistic, F = {xo ~ Xs)/iXs/-^)^ where 
is the number of degrees of freedom remaining in the data after fitting out the source. 
This statistic should be distributed according to an F-distribution with 1 and degrees 
of freedom ( [Martin 1971[ ), which allows us to calculate the cumulative probability P{F) of 



drawing a value greater than or equal to F if the null hypothesis — that the actual value 
Su{q:) = — were true. Values of F ^ 1 are highly unlikely under the null hypothesis, and 
correspond to significant detections of compact features in the data. Figure shows F{a) 
and P{F) for each of the three frequency bands. 

Since we are plotting the cumulative probability of the null hypothesis being true, 
lower values of P{F) correspond to more significant detections. We have multiple lobe 
beam patterns, so a bright unresolved feature will alias to several values of right ascension 
separated by the 40' beam throw. Thus, the clusters of peaks near 15^^ and 19'?3 do not 
correspond to multiple features, but rather the aliasing of a single feature into the side 
beams. 



We build a list of candidate unresolved features for each band by iteratively searching 
for and subtracting the most significant feature, defined here as the maximum value of 




14 16 18 20 

R.A. (J1 992.5) 



Fig. 1. — F Statistic (solid curve) and cumulative probability P{F) (dotted curve); (a) 
5.7 cm^^ band, (b) 9.0 cm~^ band, (c) 16.5 cm~^ band. Low values of P{F) correspond to 
significant detections. Note tfiat tfie multi-lobe antenna patterns cause individual sources to 
alias into several locations (see text) . 
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F{a). After each feature is subtracted, S^{a), and F{a) are re-evaluated to help suppress 
the aliasing in right ascension. We continue this procedure until -Fmax ^ H, corresponding 
to P{F) > 10~^. Three candidate features were identified this way for the 5.7 and 9.3 cm~^ 
bands. We allow each feature to be positive or negative, but the three brightest features 
are all positive. No features were found at 16.5 cm~^ which satisfy the P{F) < 10~^ 
requirement. Table |I| lists the results of this search, along with the corresponding central 
declinations. 



3.3. Spatial Extent 

To determine the spatial extent of each feature, we fit the data to a Gaussian shaped 
surface brightness model, It,e~^^~^°'^'^^'^'^^ , with the central brightness I^, right ascension 
a, and width a = ^fwhm/ (2a/2 In 2) as free parameters. The fits are performed using the 
Levenberg-Marquardt method (Press et al. 19921) . The total flux density of the feature is 



related to the central brightness by = 27ia'^I^. 

The fits near 15*^ indicate a clearly resolved feature in both the 5.7 and 9.3 cm^^ bands, 
with 6'fwhm ~ 35'. The results of these fits are given in Table ^ There are substantial 
correlations between ^fwhm and J^, so these are also listed. From I^, we can also infer an 
equivalent peak ATcmbr = 360 ± 50/iK (5.7 cm"^), 330 ± lOO/xK (9.3 cm^^), near 15^^. The 
fits near 19^'3 and 19^'9 do not show evidence for significant spatial extent. We can infer 
68% (1-0") upper bounds on 6'fwhm from the estimated parameter variance of the fit; these 
give 6'fwhm < 25' at 19^'3 and 6'fwhm < 35' at 1^9, consistent in both bands. 

Because of the close proximity of the 19^*3 and 191^9 features, there is some correlation 
between their measured fiux densities. To account for this, the fits are repeated, omitting 
^FWHM while simultaneously fitting for a and Si, of both features. The results of these fits 
are shown in Table ^. 



3.4. Limits at 16.5 cm ^ 

While there are no bright features detected at 16.5 cm^^, we can find limits for the 
corresponding features found at the lower frequencies. For the first feature, we fit using 
our extended model, fixing a = 14^^95 and ^fwhm = 35', while for the other features, we 
fit assuming completely unresolved sources with a = 19'?31 and a = 19^'92 fixed. In all 
cases, the results are consistent with zero fiux, and upper bounds can be inferred from the 
variance. These are also listed in Tables ^ and ^ 
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Table 1. Candidate Unresolved Features 



R.A.^ 


Dec* F P{F) 


S. (Jy) 




Band 1 (5.7 cm-^) 


14'?97 
19^29 
19^92 


+81?97 62 1.6 X 10-^^ 
+81?87 25 6.9 x 10"^ 
+81?85 14 1.6 X 10-^ 


+5.5 ±0.6 
+2.5 + 0.5 
+1.9 + 0.5 


— Band 2 (9.3 cm"!) — 


14'?96 
19^30 
19^94 


+81?97 27 2.5 x 10"^ 
+81?87 12 5.8 X 10-^ 
+81?85 12 4.9 X 10-^ 


+4.9 + 0.9 
+2.7 + 0.8 
+2.9 + 0.8 



^J1992.5 coordinates. 
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Table 2. Fit Results at 15^ 



Parameter Fit Result^ 


Correlations 




a ^FWHM 


5.7cin-i (X^ = 720/627) 




a : 14'?96±0^'02 




^fwhm: 36'±8' 


-0.35 


75.7: 158±24 kJy/sr 


+0.21 -0.68 


9.3cm-i (X^ = 646/627) 




a: 14'?95±0^'03 




^fwhm: 33'±13' 


-0.32 


79.3: 137±40 kJy/sr 


+0.22 -0.79 



16.5cm-i (x^ = 542/629) 

716.5: -27±77 kJy/sr 



^Errors are statistical. 
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Table 3. Fit Results at 19^3 and 19^9 



Parameter Fit Result^ Correlations 



,1^3 Q,19^9 5*19^3 



5.7cm-^ (X^ 


= 760/626) 








19^29±0^'02 








19^'92±0^03 


-0.04 




'-^5.7 ■ 


2.7±0.5 Jy 


+0.06 


-0.06 


Ql^9. 

"^5.7 • 


1.9±0.5 Jy 


+0.09 


+0.02 +0.10 



9.3cm-i (X^ = 652/626) 

a^^^^: 19^32±tf?03 

0,19^9. i9h93±o^o3 -0.12 

Slp: 3.1+0.8 Jy +0.07 -0.10 

Sl^^: 3.0+0.8 Jy +0.08 +0.06 +0.11 



16.5cm-^ (X^ = 541/628) 

Slii: 1.0+2.5 Jy 

Sli^^: 2.9+2.6 Jy +0.07 



^Errors are statistical. 
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4. Spectral Models for Features 

We now compare the measured spectra of the features with four simple models for 
astrophysical sources: CMBR anisotropy, free-free emission with 1^ oc z/~°'^, synchrotron 
emission with 1^ oc z/~°-^, and cold dust emission with oc u"^ x (i?i,(4.8K) — Bi,{Tcmbr))- 
The last of these is motivated by the possible detection of a cold dust contribution to 



the Galactic spectrum observed by the COBE/FIRAS experiment ([Wright et al. 1991| , 
Reach et al. 1995| ). Each of the models is fit to the values listed in Tables |^ and ^, the 
results of which are displayed in Figure |^. The fo^' the spectral fits are in Table ^. The 
relative calibration uncertainty is small compared to the statistical uncertainties, and is not 
included in the values. 

All three features are in agreement with the CMBR model. The first two are 
significantly discrepant with a cold dust spectrum, while the feature near 19^9 is marginally 
consistent with that spectrum. Furthermore, this location corresponds to the brightest 
warm dust emission in our field. Based on this, we cannot rule out the possibility that the 
19^'9 feature is be due to galactic dust emission. 

From these data alone, either free- free or synchrotron emission gives an acceptable fit, 
and could explain the detections. This would, however, require implausibly large amplitudes 
for emission at lower frequencies. For example, the best fit free-free spectra, extrapolated 
to 40 GHz (Q-band), imply antenna temperatures of 3.4, 1.1, and 0.9 mK for the three 
features; measurements made at Saskatoon of the same field strongly rule out such bright 
signals ([Netterfield et al. 1996|) . The limits on synchrotron emission are even stronger, 
since the extrapolated temperatures are roughly 3 times greater. Recent measurements of 
Ha emission in this part of the sky imply the free-free contribution to the MSAMl signal 
at 5.7cm~^ is <500 Jy/sr rms ( [Simonetti, 'Ibpasna, and Dennison 1996| ), several orders of 



magnitude smaller than the signal we report. 



5. Searches for Counterparts 

5.1. Catalog Searches 

We consider the possibility that these detections correspond to compact radio (or 
infrared) point sources. Note that the power law models corresponding to free-free and 
synchrotron also can be interpreted as fiat-spectrum or steep-spectrum radio spectra, 
respectively, typical of extragalactic radio sources. No counterparts for the features were 
found within one beamwidth of the feature locations in the 5 GHz S5 Polar Cap Survey 
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V (cm"^) 

Fig. 2. — Measured spectra for the three features; (a) 14^^95, (b) 19^'3, and (c) 19'?9. 
Horizontal bars show the effective width of each band. Superimposed are the best fit 
curves for a CMBR anisotropy spectrum (sohd), free- free (dashed), synchrotron emission 
(dot-dashed), and cold dust (triple-dot-dashed). 



Table 4. Spectral Fit Values 



Free- Sync Cold 
R.A. CMBR Free Dust 



151^0 


0.8 


4.5 


2.0 


40. 




0.3 


0.7 


2.0 


20. 


1^9 


2.4 


1.7 


4.3 


6.9 



Note. — All ^-re for 2 degrees of 
freedom. 
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( IKiihr et al. 198TD . For all three features, the best-fit flat spectra predict 5 GHz fluxes 
greater than 3 Jy, while the best-fit steep spectra predict fluxes greater than 40 Jy. The S5 
survey is estimated to be complete at 5 GHz down to 250 mJy, providing a strong limit on 
any flat- or steep-spectrum radio counterpart to the detections. 



A similar search was performed with the IRAS 100 /xm Faint Source Catalog (|Moshir 
et al. 19921 ). A total of 5 objects ranging in flux from 0.9-2.6 Jy were found within one 
beamwidth of any one of the three features, consistent with the average IRAS FSC source 
density of ~ 4/sq. degree. The contribution of these sources to the MS AM signal is shown 
to be negligible using the 1.25 mm to 100 /im flux ratio, (3.1 ± 0.5) x 10~^, determined 
from a complete sample of IRAS galaxies, together with a relatively flat u^'^ spectrum 
locally around 1.25 mm ( [branceschini and Andreani 1995| ). The greatest contribution is at 



16.5 cm^^, where it is estimated to be ~ 0.1 Jy for each of the latter two features, and 
~ 50 mJy for the 15*^ feature. The contributions at 9.3 cm~^ are a factor of 5 smaller than 
these, while at 5.7 cm~^ they are a factor of 30 smaller. These estimates are reinforced with 
direct observations at 90 GHz of known IRAS sources in the vicinity of the 15^ feature, 
none of which show significant flux density above the rms noise of 25-50 mJy ( Chernin and 
Scott 19951) . 



5.2. MSAMl-94 

In 1994, we flew the MSAM instrument a second time to observe the same field as the 
1992 flight. While the sky overlap between MSAMl-94 and MSAMl-92 was not perfect, it 
was sufficient to permit an overall comparison of the measurements ( [Inman et al. 199B ). We 



have repeated the fits of §^]3|on the MSAMl-94 data to see if the bright features are present. 
Near 15*^, there is a corresponding hot spot at a = 14'?85 ± 0'?03 with ^fwhm = 38' ± 12', 
Is. 7 = 106 ± 32 kJy/sr, and Iq ^ = 71 ± 62 kJy/sr (statistical errors only). Near 19^*3 and 
19'^9, there are corresponding features that are fit simultaneously as before. At 5.7 cm , 
the resulting parameters are a = 19^31 ± 0'?02, Slf^ = 1.3 ± 0.3 Jy, and a = 19^93 ± 0'?02, 
^19^9 = 1.4 ± 0.4 Jy. At 9.3 cm'^, the fit gives a = 1^26 ± 0^03, S^^^ = 1.8 ± 0.5 Jy, 
and a = 19'?93 ± 0'?02, S^p = 3.2 ± 0.7 Jy. The correlations in the MSAMl-94 fits are 
larger than the MSAMl-92 fits; in particular, the 1^3 and 1^9 flux densities have variance 
correlations of +0.49 at 5.7 cm~^, and +0.42 at 9.3 cm~^. The somewhat smaller flux 
densities in MSAMl-94 are consistent with the pointing differences between the two flights. 
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5.3. Other Observations 

Since the initial report of these detections in Paper I, a number of groups have searched 
one or more of these regions for corresponding sources at different wavelengths and/or 
resolutions. Of particular note is the search reported by [Church et al. 1995| using the 
Caltech Submillimeter Observatory on Mauna Kea. This search, which made a deep map of 
the sky around the 15^^ feature at 4.7cm~^ using a V.7 beam, places an upper limit on any 
point source in that field of 5*4.7 < 1 Jy- 

Another observation of this field has recently been reported by [Netterheld et at. 1996 . 
These observations were made from the ground at Saskatoon, Saskatchewan, Canada, in 
Q-band (36-46 GHz, 1.2-1.5 cm^^) at 0?5 resolution, a portion of which were designed 
to reproduce, as closely as possible, the MSAM observations. A three-lobe beammap, 
similar to the MSAM double-difference demodulation, was synthesized (the single-difference 
demodulation is not usable from the ground due to atmospheric noise). Although having 
somewhat less statistical weight than the MSAM observations, the Saskatoon results 
correlate well with MSAM, including the bright feature near 15^ right ascension. 



6. Simulations 

The presence of bright, compact features in the MSAMl-92 data was unexpected. 
In Paper I, we questioned whether CMBR fluctuations obeying Gaussian statistics could 
produce such features. To address this issue, we have performed a detailed Monte Carlo 
simulation of the MSAM observations. 

The approach used is similar to that outlined in §4.3 of Paper iQ. We generate numerical 
realizations of the observations, including contributions from instrument noise and CMBR 
temperature anisotropics, assuming the anisotropy is a Gaussian random field described by 
a two-point correlation function, C{9). In Paper I, we set limits on the overall amplitude 
of l\T jT = \JC{0) by using the likelihood ratio statistic A; here, we use those upper and 
lower bounds to generate realizations which are then searched according to the procedure 
of §3.2| for unresolved features. Unlike Paper I, here we simulate both demodulations 
simultaneously, using 6c = 0?4. We thus determine the distribution of n, the number of 
unresolved features with P{F) < 10^^, corresponding to the range of AT/T fluctuations 
found in the data. Our measured number of sources is n = 2 or 3, depending on the 



^Eq (2) in Paper I contains a typographical error, and should read 

(sfcS;) ^ J dXidX2Hk{xi)Hi{x2)C{\xi - X2\). 
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assignment of the 19^'9 feature to CMBR or galactic dust. For AT /T = 1 x 10 ^, we find 
the cumulative probability P{n > 2) = 0.20, P{n > 3) = 0.04, while for AT/T = 3 x 10-^ 
P{n > 2) = 0.90, P{n > 3) = 0.69. From this, we conclude that bright, compact features 
such as those observed are in fact consistent with our assumption of Gaussian statistics. 

We have also studied the distribution of n using a more physically motivated 
correlation function, calculated for "standard" flat cold dark matter (CDM) with 
Vt = 1,Qb = 0.05, A = 0,Hq = 50km/s/Mpc, and no reionization ( |Sugiyama 1995| ). Using 



this model, normalized to the COBE/DMR detection Qrms-ps = 20/iK ( porski et al. 199^ ), 
we find reasonable probabilities for the observed n: P{n > 2) = 0.43, P{n > 3) = 0.15. 
Thus, the detections are consistent with adiabatic fluctuation models such as CDM. Similar 
results have also been reported for a more generic study of 0?5 resolution CMBR maps 
( [Kogut, Hinshaw, and Bennett 1995| ). 



7. Conclusions 

All three features are best described as temperature anisotropics in the CMBR, 
although the possibility remains that one of these is due to galactic foreground dust 
emission. While we were initially surprised to find such bright, compact features in the data, 
we now understand that these are consistent with our assumption of Gaussian fluctuations. 
Although our catalog searches were all null, we find the most compelling argument that 



the features are not true point-sources comes from the direct search by [Church et al. 1995 
Finally, the morphological agreement reported by the Saskatoon experiment, consistent 
with a CMBR anisotropy spectrum spanning three octaves in frequency from Q-band (1.2 
cm~^) to the MSAM 9.3 cm~^ band, leaves little room for alternate explanations of the 
detection. Given the agreement with a CMBR temperature anisotropy spectrum, together 
with the implausability of other explanations, we conclude that at least the first two, and 
possibly all three of these detections, are cosmological in origin. 
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